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1. INTRODUCTION 

Congestion is an important issue in power system operation. Transmission congestion occurs when a 
power transmission system is unable to meet the needs of consumers due to higher demand. This gives an 
adverse effect on the electrical energy market where it cannot operate at a competitive level and may 
endanger the safety of the system [1]. Congestion management is a solution to address the issues of the 
operating system safely and consistently in an electric power system. The incidence of congestion can occur 
in several forms. Among these are caused by the failure of the transmission line, interference generator, 
changes in demand and clumsy dealing. The impact of congestion in a power system can lead to increased 
costs in certain areas and damage the device in the system [2]. 

Voltage stability phenomenon is a progressing issue in power system. It is one of the subsequent 
events which may occur due to congestion. It is normally measured using a stability index. C. Subramani 
in [3] proposed a Global Voltage Stability Index for voltage stability assessment. This index was employed 
with Artificial Neural Network to solve the voltage collapse problem. Since its introduction, ANN has 
become very popular, and applied to various decision making problem in numerous topics such as fault 
location [4-5], frequency, online ATC and real time power system estimation [6-8] voltage security 
evaluation [9], prediction of geothermal power plant system performance [10], and the rest of researchers 
was discussed about voltage stability monitoring [4, 11-13]. Today’s scenario has faced of steady increase in 
demand of electricity. With the development of different FACTS devices, compensation effort has become 
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another option for power system operators. Hence there is a need to study the optimum placement of FACTS 
devices. This is to avoid the over-compensation or under-compensation phenomena if no optimization 
process is involved. Various models for SVC, TCSC and PSS are conceivable and applied in different studies 
[14-22]. According to this issue, K. Vetrivel et al [23] suggests the implementation of Thyristor Controlled 
Series Capacitor in order to solve the congestion problem using Genetic Algorithm (GA) technique. In [24], 
the researcher highlighted a fuzzy technique to determine the optimal location of Thyristor Controlled Series 
Capacitor (TCSC). Therefore, the researcher used a Line Utilization Factor (LUF) and Real Power 
Performance Index (RPPI) factor to find out the level of congestion in a transmission line. Lastly, to observe 
the presence of TCSC, sensitivity parameters for the total system losses are derived as a function of the real 
power at the individual load points. 

In this paper, Integrated Multi-Layer Artificial Neural Networks (IMLANNs) model has been 
developed for congested line prediction in choosing power system. Then, a new optimization technique 
termed as Clonal Evolutionary Particle Swarm Optimization (CEPSO) was developed in order to determine 
the optimal sizing of single unit of SVC and TCSC. The results highlight the advantage of new technique as 
compared to other technique with the objective function to minimize the transmission loss reduction and 
improvement of voltage profile in the system. 


2. RESEARCH METHOD 
2.1. Integrated Multi-Layer Artificial Neural Networks (IMLANNSs) 

Firstly, this research started with involves the development of a new model for determining the 
location of a FACTS device using artificial intelligence technique: namely Integrated Multi-Layer Artificial 
Neural Networks (IMLANNs) technique for a single objective function for predicting the Fast Voltage 
Stability Index (FVS/) value and congested lines in power system. For the purpose of validation, the proposed 
technique is tested on two test systems, namely IEEE 30-Bus Reliability Test Systems (RTS). Before the 
IMLANNs model was developed, a case study was chosen to get an early result in power system load flow 
during normal and contingency condition. A pre-developed index, FVSI/ is used as a tool to identify the 
congested lines [25]. The experimental results of FVS/ value is chosen as the targeted output in the 
IMLANNs model. 





2.1.1. Algorithm for the Proposed IMLANNs 

This section presents the whole procedures of the proposed IMLANNs model. The process of the 
proposed IMLANNs is explained below: 
1. Choose the load bus. In IEEE 30-Bus RTS, there are 11 load buses that have been chosen for the testing 
process. 
Read system data: IEEE 30-Bus RTS. 
Set the loading condition. The system allows to run at the maximum permissible load. 
Run load flow. 
Calculate FVS/ values for all lines and identify the congested lines. The equation of FVSI is show in 
(1). 


oe 


42? Or 
FVSI = ee () 





6. Build the raw data patterns for the ANN. If the number of patterns is not adequate, adjust the loading 
condition. 

7. Cluster the training patterns and testing patterns. 

Call the training data patterns. 

9. Initialize the training parameters. On the stage of the training process, the numbers of nodes in each 
hidden layer, momentum rate, learning rate and other training parameters have been conducted by 
extensive experiments in order to determine the optimal value. 

10. Develop the training configuration. The number of neurons in the first hidden layer, the number of 
neurons in the second hidden layer, the type of active of activation functions, number of iterative 
updates (epochs) and the mean square error goal can be varied in the process of determination optimal 
types or values. 

11. Perform the training data patterns. 

12. If the solution does not converge, vary the training parameters. 

13. Save the IMLANNs model. 
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14. Call the testing data patterns. Retrieve the saved IMLANNs model. The testing process is conducted 
consecutively after obtaining the best trained ANN model. 
15. Perform testing process. Calculate the absolute error (E,»;) as follows: 


Fabs = |o, a 0p, | (2) 


16. Where O, is the calculated output and QO, is the predicted output. n is the number of training patterns 
involved. The regression analysis is used to minimize the prediction error when fitting the predicted 
output to the targeted output. 

17. If Ris not close to 1.00, repeat step xii. 

18. Tabulate the IMLANNs output and targeted output. 


2.2. Clonal Evolutionary Particle Swarm Optimization (CEPSO) 

Consequently, a hybrid optimization technique, namely Clonal Evolutionary Particle Swarm 
Optimization (CEPSO) is developed in order to determine the optimal sizing of single and multi-unit of 
FACTS device installation. The idea of hybrid CEPSO is developed based on the combination between 
Artificial Immune System (AIS), Evolutionary Programming (EP) and Particle Swarm Optimization (PSO) to 
improve the weakness experienced in the conventional PSO. This technique was built up to solve single 
objective problem (SOP) which considers minimizing the power loss and to improve the voltage profile. 


2.2.1. Algorithm of CEPSO Technique 

In this study, CEPSO is implemented in order to minimize the total losses in the distribution system; 
while monitoring voltage at all buses in the system to be within the acceptable limit. The step by step CEPSO 
algorithm for the proposed optimal sizing of FACTS devices is explained as follows: 
1. Set the PSO parameter (Table 1). 


Table 1. PSO Parameter 








Paramater Value 
Total number of particles, Np 1,2,3 
Maximum iteration, Kmax 100 
Minimum weight, Winin 0.4 
Maximum weight, Wimax 0.9 
Learning rate for individual ability, C; 2.05 
Social influence, C2 2.05 
Velocity constant factor, € 1 
Minimum particle position, Ximin 0 
Maximum particle position, Ximax 1 





Set the parameters of the control variable to optimize as sizing of SVC and TCSC. 

Set the loading condition for pre-optimization. Set the loss and voltage constraints, i.e. preLoss > 
postLoss and preVoltage < postVoltage, where preLoss and preVoltage are the loss and voltage profile 
values obtained by running load flow before any compensation process is conducted. This is important 
to certify that all generated initial population satisfy all the equality and inequality constraints. 

4. Generate an initial population to consider the variable that should be optimized; the sizing of single or 
multi FACTS device. The random numbers x is the control variables for the FACTS device (x;, 
Xp....X,) Where x, is the sizing of FACTS devices to be installed into the system. The general matrix for 
initial population; X,,,, and X¢¢s¢ is given by (3) and (4) respectively. 


Py 


Xsve = [%1 x2 Xn] (3) 
Xtesc 7 [%1 x2 Xn] (4) 


5. Initialize velocity and position. 

6. Calculate Fitness 1. By considering the objective function, fitness is the particle in the population, 
which is evaluated and calculated by running the load flow. 

7. The Clone process is a process to duplicate the parents. The mathematical equation is given in (5). 
Cloned Population, Timk:- 
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My ee Mae ot Xa kt. X1kt20  %1,2k l 
: : i 
Xn1 oo Xnk os Xntak Xak+20°  Xn,2k 
Xo1 se X2R oe XO kt. Xak+2e+ %2,2k 
Tamk = : m (5) 
Xnt oo Xnk oo Xnkt1 Xak+20  Xn,2k 
X31 0 Xn2  AXnkt+1 Xaet2- *Xn,2k l 
7 : m 
Xn1 oss Anz + Ankt1 *X1k+20  Xn2k 


Matrix size of cloned population: mn x 2k 

where: 

n = population number 

m = cloning number 

k =number of control variables 

8. Calculate the new particles of position using mutation technique given in (6) and (7). Mutation is a 
process to produce offspring (children). £ is a scalar value, generated randomly to control the mutation 
process. 6 controls the movement of candidates from the valley to hill within the chosen bell-curve. B 
can be pre-determined or randomized. For randomized B value, this technique will become adaptive. 
This curve can be the Gaussian, Cauchy or Levy. The matrix sizing for the offspring is similar to that 
for the cloned population, Tpmx- 





Xi+m,j = Xij + N(0, y’) (6) 
* = BC maxX min) (2 

you p Xjmax Xjmin fax ( ) 

where: 

0<B <1 =rand (1,1)*1 

Bb = mutation scale 

fi = fitness for i” random number 

Triad = maximum fitness 

N = Gaussian random variable 

Xij = parents 

Xfmax = maximum random number for every variable 

Xpmin = minimum random number for every variable 

Xivmj = mutated parents (offspring) 


9. The combination is a process to group together the whole population and population after the cloned 
process in cascode form. It is can be conceptually represented as in (8), (9) and (10). From the parent 
population, A,and offspring population, Az. Therefore, the combined population can be written as C. 














NOX, XQ XK Xt Xk+2""* *2k Fy 
Aj = : st as : : ras : : (8) 
n Wee. Waewes wee ree Wee ee Xn2kFo9 
no Xinew *X2new'’’ %iknew *1k+new eh Fy 
A, = : : og : : io: (9) 
nxm eee eee see eee eee °F oq 
Ay 
c=|4| 10 
a (10) 


10. Tournament selection is a process to prescribe the candidates for the next iteration. If the cloning 
multiplier m is 10, therefore Fl and F2 will have 20x10=200 populations. Only 20 best members or 
individual is prescribed from this population for the next iteration. There are many techniques for the 
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selection process such as pairwise comparison, elitism and a roulette wheel. Any suitable technique can 
be adopted for this purpose. But in this study, a pairwise comparison was used. 

11. Find the Prest_ oid ANd Gpest ola Value and stored as Prest oig= MIN (Phest fitness) old ANA Goest_old= 
Prest_fitness_min_old. 

12. Update the velocity and position of the particle.The condition to update velocity and position are 
described as follows: 


Update velocity, If VES Vine 
set vis = Vnax 

Else if V;* < Vinin 
set VK = Vinin 


Update position, If RES Res 
Set Ke = aos 

Else if x* < Xmin 
Set = xi 


13. Calculate Fitness 2 and determine the Prest new ANd Grest_new Value. Prest new = MIN (Phest_ fitness) new and 
Ghest_new = Prest Fitness_min_new. 

14. Perform convergence test. The convergence criterion is defined as the difference between Fitness and 
Fitness 2. Repeat steps (vil-ix) until stopping criterion, as such sufficiently excellent Gye fitness or a 
maximum number of iteration is met. 

15. End the CEPSO process. 


3. RESULTS AND ANALYSIS 

Firstly, loss minimization was chosen as the first objective function in this section. All the SVC 
locations were determined from the IMLANNs technique and the optimal sizing of single unit of SVC 
installation using the CEPSO technique in the IEEE 30-Bus RTS are tabulated in Table 2. Five load buses 
were chosen with their maximum reactive power loading in order to observe the installation of SVC to 
manage the congestion in power system. The participating load buses are buses 3, 4, 7, 10 and 12. With 
maximum reactive loading, the congested line was identified based on the calculation of FVS/ as the 
benchmark of congestion situation. From the congested line, the bus connected from bus to another bus can 
be known. The location of SVC was installed at the bus which experiencing congestion. 


Table 2. Results of the single unit SVC installation using the CEPSO technique in IEEE 30- Bus RTS 








Loaded Qumax Optimal - Loss 
bis (MVar) SVC Loc SVC Sizing Pre-opt Post-opt (%) A Loss 
(Bus) (MVar) (MW) (MW) 
3 351.2 4 10 69.7726 65.1713 6.59 
4 391.6 12 10 49.4885 48.9081 1.17 
7 350.9 7 10 83.3981 76.9638 7.72 
10 182.0 11 10 30.1015 28.3162 5.93 
12 187.5 12 10 21.9468 21.5076 2.00 





(%) A Loss= % of loss reduction 


At Q47=350.9MVar, Bus 7 has been determined as the SVC location while the optimal SVC sizing is 
10MVar identified by using CEPSO. The percentage of minimum loss is 7.72%, which is significantly high, 
compared to other load buses. For instance, at load Bus 3, the percentage of loss minimization show the 
second higher worth 6.59%. The optimal SVC sizing is 1OMVar while Bus 4 has been determined as the 
SVC location. From the observation, the overall result was shown the excellent reduction of loss 
minimization when the SVC was installed as a single unit in IEEE 30-Bus RTS. 

TCSC was chosen as the second FACTS devices to manage congestion in power system. TCSC was 
installed on the line which apparently identified as congested line. The optimal location and sizing of TCSC 
was identified using CEPSO. Table 3 presents the results from the single unit TCSC installation using the 
CEPSO technique on the IEEE 30-Bus RTS. The table tabulates the optimal location and sizing of TCSC. 
The losses before and after CEPSO was implemented on the system are also shown. The percentage loss 
denotes the % of loss reduction with CEPSO implementation. On the overall, the power loss with TCSC 
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installation has significantly reduced the loss for each loaded bus as shown in the table. The highest 
percentage of loss reduction of 3.14% when Line 15 was identified as the suitable TCSC location with its 
corresponding optimal sizing of 1.5531p.u.. At Qyjo=182.0MVar, the optimal TCSC sizing is 0.6897 p.u., 
while the TCSC location is Line 13. At this condition, the percentage of loss reduction is 2.29%. For optimal 
TCSC sizing at Qai2 of 187.5MVar, is 0.3948 p.u. with the optimal location is Line 16. The percentage of 
loss reduction is 1.79%, reduction from 21.9468 MW to 21.4118 MW. 


Table 3. Results of single unit TCSC installation using CEPSO technique in IEEE 30-Bus RTS 








Optimal Loss 
mee aa TCSC Loc TCSC Sizing Pre-opt Post-opt oe 
(Line) (p.u.) (MW) (MW) 
3 351.2 4 0.3674 69.7726 69.4052 0.53 
4 391.6 15 1.5531 49.4885 47.9354 3.14 
7 350.9 8 1.4213 83.3981 81.9768 1.70 
10 182.0 13 0.6897 30.1015 29.4118 2.29 
12 187.5 16 0.3948 21.9468 21.5520 1.79 





(%) A Loss= % of loss reduction 


In order to monitor the voltage profile improvement, the results of single unit SVC installation in 
IEEE 30-Bus RTS are illustrated as bar graph in Figure 1. In this figure, observation can be seen that the post 
voltage shows the voltage improvement as compared to the pre-voltage to all location of single unit SVC in 
the system. As mentioned in previous section, five SVC locations are chosen from the loaded buses. These 
locations are buses 4, 12, 7, 11 and 12 as seen on x-axis. For example, voltage after single unit of SVC is 
installed at Bus 4 has been improved to 1.0236p.u from its original value of 0.8264p.u. The detail results of 
all single unit SVC installation are depicted in the figure. 


Single Unit SVC 
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Figure |. Results of voltage improvement when single unit SVC installation in IEEE 30-Bus RTS 


Next study is the installation of the single unit of TCSC in the IEEE 30-Bus RTS. The results of 
voltage profile improvement when single unit of TCSC installation in IEEE 30-Bus RTS are presented in 
Figure 2. In this figure, observation can be seen that the post voltage shows the improvement as compared to 
the pre-voltage to all location of single unit TCSC in the system. As mentioned in the previous section, five 
TCSC locations are chosen from the loaded buses. These optimal locations are lines 4, 15, 8, 13 and 16 as 
seen on the x-axis. For example, voltage after single unit of TCSC is installed at Line 13 has improved to 
0.9893p.u. from its original value of 0.9591p.u.. On the other hand, voltage before installation the single unit 
TCSC is increased when Line 15 is chosen as location of TCSC from 0.8632p.u. to 0.941 1p.u.. The detail 
results of all single unit TCSC installation are depicted in the Figure 2. 
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Figure 2. Results of voltage improvement when single unit TCSC installation in IEEE 30-Bus RTS 


4. CONCLUSION 

This paper has presented a new two technique as IMLANNs and CEPSO technique. IMLANNs 
technique was developed to predict the congested line in the system while CEPSO used to solve the single 
unit installation of SVC and TCSC for two independent objective functions. The first one is to minimize the 
loss; while the second objective is to improve the voltage profile. Both objectives have been successfully 
achieved through the implementation on IEEE 30-Bus RTS. The important effects to be highlighted in the 
results is the minimization of loss with the single unit of SVC and TCSC installation into the system. From 
the results, it is observed that the installation of single unit of SVC and TCSC with the optimal locations and 
sizing can help the power system operator for their future planning. 
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